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Abstract 



I analyse the model of Babu and Mohapatra for the fermion mass matrices in minimal 
SO(IO). Those authors have not considered the whole variety of fits to the experimental 
data possible with their mass matrices. Consequently, their predictions for the neutrino 
Q> . spectrum are incomplete. I survey various types of neutrino spectra possible in their 

fNl ! model. 

O 

m 

^ I Babu and Mohapatra (BM) |^ have recently noted that in a minimal SO(IO) model the 

^ ' lepton mass matrices can be written as functions of the quark mass matrices: 

iu: Ml = aMu + bMd, 

= M-2)M„ + ^^±^^1^M., (1) 



Mu and Md are the mass matrices of the up-type and down-type quarks, respectively. Mi is the 
mass matrix of the charged leptons. and M^ are the Dirac and Majorana mass matrices, 
respectively, of the neutrinos. The effective mass matrix of the light neutrinos is, in the seesaw 
approximation, 

Mf = -Ml^M^\Ml^f . (2) 

All the mass matrices, in Eqs. ^ as well as in Eq. ^, are real and symmetric, a and b are real 
numbers.^ 

The model of BM is attractive because of its minimality in the content of scalar represen- 
tations. BM only use one 10 and one 126 of Higgs scalars coupling to the fermions. These 
representations are necessary anyway: the 126 is needed to give Majorana masses to the neu- 
trinos and to generate the seesaw mechanism; the 10 is needed in order to obtain nib ~ at 

^The scheme proposed by BM is more general, in that it also includes two phases, which render the mass 
matrices complex. But BM have only analysed in detail the case of real mass matrices and I shall also consider 
only tliat restriction of tlieir scheme. 
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the breaking scale of the Pati-Salam subgroup SU(4)p5' of SO(IO). The predictive power 
of the BM model originates in the small number of Yukawa-coupling matrices, not in ad hoc 
assumptions, like discrete symmetries leading to vanishing Yukawa couplings, which are fre- 
quent in other models BM have noted that the 126, besides giving Majorana masses to 
the right-handed neutrinos, also introduces welcome corrections to the other mass matrices. 
Therefore the BM model is worth especial attention. My purpose in this Brief Report is to fit 
the known data on the quark masses and mixings, and on the charged-lepton masses, by mass 
matrices of the BM type, and then to find out which kinds of neutrino spectra are possible in 
the model. 

The BM model leads to two conditions on the quark masses and mixings, and on the 
charged-lepton masses. The first condition is that, as the quark mass matrices Mu and Md are 
real, the quark mixing matrix V is real. It is well known p that, once \Vus\, \Vcb\, \Vub\ and 
\Vtd\ are given, the quark mixing matrix is completely determined. I have used the following 
fixed values in my fits: \Vus\ = 0.22, \Vcb\ = 0.046, and \Vub\ = 0.004. The condition that 
V is orthogonal then leads to two possible values for \Vtd\- either \Vtd\ = Knin = 0.006223 or 
l^dl = Knax = 0.014017. (For Vmin < \Vtd\ < Knax the matrix V is unitary but not real; for 
\Vtd\ < Knin or \Vtd\ > Knax the matrix V is not unitary.) If one accepts the standard- model 
explanation of Bd-Bd mixing by means of the box diagram [Q and uses a top-quark mass of 
150 GeV (see below), then \Vtd\ = Knax corresponds to fs = fB\/BBVB = 116 MeV, a quite 
reasonable value. \Vtd\ = V^in corresponds to fs = 262 MeV, a rather high value, but still 
acceptable. 

The second condition can be derived from the first of Eqs. Obviously, 
trM; = a tiMu + h trM^ , 

trMf = iiMl + 2ah tr (M^M^) + 6^ trM^ , (3) 
trMf = iiMl + 3a^6 tr(M^Md) + 3ab'^ tr(M„Mj) + trMj . 

Because M„ and Md are Hermitian, tr(M^Mj') = J2i=u,c,t Y.j=d,s,bf^^if^]\^ij\'^ ■ Therefore, all 
the traces in the right-hand sides of Eqs. ^ are observables. The traces in their left-hand sides 
are functions of the charged-lepton masses. We thus have three equations for the two unknowns 
a and h. From the first equation, h is determined as a function of a. Substituting the result in 
the second and third equations, we obtain a set of a quadratic equation and a cubic equation 
for a: 

xio^ + 2x2a + X3 = , 
l/ia^ + 32/20^ + 3?/3a + 2/4 = 0, (4) 

where the coefficients Xi^2,3 and 1/1,2,3,4 are functions of traces, and are therefore observables. 
Eqs. H fix a to be 

3x1X31/2 - 2x2X:>,yi - xly^ 

= TT^ 7. • yp) 

Axivi + 3a;f?/3 - QxiX2y2 - xiXsyi 

Besides, because one has two equations for only one unknown, the coefficients of Eqs. ^ must 
satisfy the following condition: 

xlyj + 9xlxsyi - 6x^x21/31/4 
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-6XiX3l/2Z/4 + 6X1X2X31/17/4 - 1 8X1X2X37/21/3 
+ 9X1X37/2 - 6X1X37/17/3 + 12XiX27/2l/4 

-6x2X37/17/2 - 8x37/17/4 + 2:37/1 + 12x2X37/17/3 = . (6) 

This is the second condition on the quark masses and mixings, and on the charged-lepton 
masses, to which I referred above. 

I have used fixed and exact values for the following masses: rrit = 150 x 0.506 GeV, nib = 
4.25 X 0.327 GeV, = ±1.27 x 0.286 GeV, = ±5.1 x 0.273 MeV, = 1.784 x 0.960 
GeV, m^, = ±105.6584 x 0.960 MeV and = ±510.999 x 0.960 keV, in which I have used the 
renormalization factors given by BM. The masses of the top and bottom quarks may, without 
loss of generality, be set positive. The masses of all other quarks and of the charged leptons 
may have either sign. However, I am interested in a picture of the mass matrices in which the 
contributions from the 10 of scalars dominate over the contributions from the 126. For every fit 
of the data, one may make the transformation a —>■ —a and 6 — > — 6 to obtain another fit, which 
has symmetric charged-lepton masses, but a different neutrino spectrum and a different leptonic 
mixing matrix. For each particular fit, one may separate the contributions of the 10 and of the 
126 to Md and Mi: the contribution of the 10 is {3Md + M/)/4, and the contribution of the 
126 is {Md — Ml) /A. I have computed the ratio of the highest eigenvalue of the contribution of 
the 10 to the highest eigenvalue of the contribution of the 126. I have found that, if mb and 
rrir have opposite signs that ratio is close to one, and I have discarded those fits. But, if and 
ttIt- have the same sign, the ratio is typically somewhere between 5 and 15; those are the fits 
that satisfy my (and BM's) prejudices. Thus, from each pair of fits related to each other by the 
symmetrization Mi —Mi referred to above, I only consider the fit in which irir is positive. 

My fitting method has been the following. Having fixed all the quark masses, together with 
their signs, except ms and tti^, and having I looked for pairs of values of ms and md such that 
the condition of Eq. ^ is satisfied. For each such pair of values, one automatically obtains a 
fit of the experimental data by the BM mass matrices, i.e., values of a and h. I was careful to 
restrict my search of and to the one-standard-deviation ranges suggested by Gasser and 
Leutwyler 0: 18.0 < \ms/md\ < 21.2 and 1.63 < \md/mu\ < 1.89 at the renormalization scale 
1 GeV. (Remember that |m„(l GeV)| = 5.1 MeV is kept fixed.) I also used the renormalization 
factors for and rrid suggested by BM. 

I found four different types of fits. All four types of fits have a positive muon mass and a 
positive ratio nis/md. Fits of type 1 have positive charm mass, all others have rric < 0. Fits of 
type 2 have positive strange mass, all others have Trig < 0. Fits of type 4 have \Vtd\ = Knax, all 
others have |Vfd| = Knin- For all four types of fits, the signs of the up mass and of the electron 
mass are not very relevant, solutions being found for any choice of those signs, and the resulting 
neutrino spectra being very similar. For each type of fit and for each choice of signs for rrie and 
rUu, a curve may be drawn in the mg vs. rud plane, displaying the pairs of values for which the 
fits are found (which yield solutions of Eq. ^). The fits presented by BM in their original work 
are all of type 4. 

I emphasize that I have been very exacting in my assumptions. I have not allowed any 
quark masses, except for and rrid, to vary, and I have used the central values suggested by 
Gasser and Leutwyler; I have not allowed the renormalization factors to vary from the values 
suggested by BM; and I have considered only the one-standard-deviation ranges for Im^/mrfl 
and \md/mu\. I have checked that if one allows instead those two ratios to be anywhere in their 
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three-standard-deviation ranges, one obtains other types of fits to the data, beyond types 1 to 
4, but I do not consider those other types of fits in this Brief Report. I also fixed V except for 
the ambiguity \Vtd\ = Vmax or \Vtd\ = Vmm, and I have discarded the fits with {mbmr) < 0. 

In the fits of type 4 the variable r2 introduced by BM may be infinite (which corresponds 
in my notation to 6 = —3), indicating a null contribution of the 126 to Mu and to . Such 
a possibility does not arise in the fits of other types. 

I now turn to the neutrino spectra. Having found a and b one may compute M^^ . Its 
eigenvalues m^, 1712 and mi, ordered as Im^l > \m2\ > |mi|, are proportional to the light- 
neutrino masses. From the diagonalization of Mi and of Mf^ one computes the lepton mixing 
matrix K. I find that the four kinds of fits yield four different types of neutrino spectra. 

Fits of types 1 and 2 yield well-defined neutrino spectra, with small variations of the various 
parameters as one lets and evolve along the curves giving the solutions to Eq. K is 
dominated by a small first-second generation mixing in both cases, and the ratios of neutrino 
masses are small. For fits of type 1 we have 



23 < Img/mal < 27 
1.6 < |m2/mi| < 2.1 
0.325 < \K^,f,\ < 0.385 
0.045 < \K^,^r\ < 0.054 
0.031 < \K^^r\ < 0.037 



(7) 



For fits of type 2 we have 



17.9 < Ima/msl < 18.8 
11.2 < Ims/mil < 12.7 
0.091 < \K^,^\ < 0.104 
0.010 < \K^^r\ < 0.012 
0.014 < \K^^r\ < 0.016 



(8) 



In contrast with these, fits of types 3 and 4 lead to less well-defined neutrino spectra. Fits 
of type 3 are characterized by very large, even maximal, lepton mixing. They yield 



1 < \m^/m2\ < 12 
24 < Ims/mil < 33 
0.35 < \K^,^\ < 0.75 
0.09 < \K,r\ < 0.13 
0.24 < \K^^r\ < 0.71 



(9) 



When one lets and evolve along the curves giving the solutions of Eq. ^ the mixing 
of the first two generations, given by If^^j^l, diminishes, while simultaneously the mixing of 
the last two generations, given by \K,y^r\, increases. Also, \m2/mi\ increases while \m^/m2\ 
decreases, and this in such a way that their product remains approximately constant (see the 
second of Eqs. ^). At a certain point the masses of the two heaviest neutrinos become equal, 
at which point the second and third rows of K must be interchanged. 
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Fits of type 4 were the only ones considered by BM. They are characterized by large ratios 
among the neutrino masses. One has 



700 < Img/mal < 2600 
20 < Ima/mil < 30000 
< \K,,^\ < 0.4 
0.006 < \K^^r\ < 0.060 
0.103 < \K^^^\ < 0.111 



(10) 



Let us confront the various neutrino spectra with the MSW explanation of the solar- 
neutrino deficit. Starting by the solutions of type 4 (see Eq. [1^), we notice that the large 
second-third generation mixing allows us to make use of the experimental data on v^-Vj. mixing 



rO| to obtain \m^\ < 2 eV. Then, |m3/m2| must be smaller than 1300 if \m2\ is not to be 



too small for the MSW effect to explain the data. This is possible for solutions of type 4, and 
therefore it is possible to accomodate a Ug-i'^ MSW effect in some of these cases [|[|. But, if this 
occurs, future experiments (P860) should be able to detect z/^-i^r oscillations, because mj,^ must 
be very close to 2 eV. Another viable possibility is z/g-z/,- MSW depletion of the solar-neutrino 
flux, in which case however all the neutrino masses would be very small, and uninteresting from 
the cosmological point of view. If MSW Ue^^^r oscillations explain the solar-neutrino depletion 
in this type of solutions, then the z/^ mass will be in the range 10^^^-10~^^ eV, and this would, 
by means of normal vacuum oscillations, provoke a further depletion of the neutrinos observed 



in the gallium experiments |11 



The MSW effect cannot occur in solutions of type 2 (see Eq. |), because the i^i-us mixing 
is too small while the mixing falls in the range between the small-mixing-angle and the 

large-mixing-angle fits of the MSW effect to the data. Solution 2 is interesting because of its 
very small i'2^^3 mixing, which is such that Im^l is unconstrained by the existing results on 
VuT^^T oscillations. The situation will change when the future results of CHORUS, NOMAD 



and P803 |T0| become available, which results will be able to exclude or confirm ranges of 
Imsl interesting from the cosmological point of view, like \m^\ ~ 5 eV (hot dark matter) or 
\m^\ ~ 50 eV (closure of the universe). Meanwhile, \m^\ can be indirectly constrained: because 
sin^(26'e/^) ^ 0.04, the E776 results ||12| imply rriy^ < 0.6 eV, and as we know the ratio |m3/m2| 
to be around 18, we obtain m^^ < 11 eV. This already excludes a substantial contribution of 
the light neutrinos to the closure of the universe. 

Solution 3 (see Eq. H) easily accomodates the large-mixing-angle MSW solution. In spite of 
the very large first-second and second-third generation mixings, a three-generation treatment 
of the MSW effect [13| does not seem essential, because l-ft'iyirl remains small. It is also possible 



that m^^ 10"^ eV leads to MSW suppression of the GALLEX and KAMIOKANDE neutrinos 
while rriy ^ 10~^ eV suppresses the Homestake neutrinos, with maximal z/g-t'^j mixing [|l4 



The large second-third generation mixing implies that \m^\ must be smaller than 0.7 eV, with 
prospects of tightening this bound to 0.25 eV in the P860 experiment [1^. This reasoning is 
invalidated if Im^l is sufficiently close to \m2\: in some solutions of type 3 those two neutrino 
masses happen to be exactly equal. But then the large mixing of z/g with z/2 and z/3 would imply 
anyway, via the Goesgen results on missing z7^ PI, that Imsl < 0.2 eV. 



Solutions of type 1 (see Eq. |^) can sometimes accomodate the large-mixing-angle MSW 
solution, in which case one obtains |m3| pa 0.13 eV, |m2| ~ 0.005 eV and |mi| ^ 0.003 eV. 
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An alternative possibility is the small- mixing- angle MSW effect between v,. and v^-. Then, 
m| «i 4 X 10~^eV^ and «i 6 x 10~^eV^, and the v^-v^^ MSW effect induces an important 



additional suppression of the neutrinos observed in the gallium experiments [14]. In any case 



the neutrino masses will be cosmologically very small, and untestable by means of terrestrial 
oscillation experiments. 

In conclusion, the BM model is not as predictive as its authors have suggested, because 
various fits to the known data are possible. Various different possibilities for the neutrino 
masses and for the lepton mixing arise. Queer possibilities, such as very small neutrino mass 
ratios or very large lepton mixings, may occur. 

I acknowledge discussions with Professor L. Wolfenstein, who read the manuscript and 
made many useful suggestions on the analysis of the neutrino spectra. I thank Professor R. N. 
Mohapatra for a short e-mail discussion on this work. This work was supported by the United 
States Department of Energy, under the contract DE-FG02-91ER-40682. 



References 

[1] K. S. Babu and R. N. Mohapatra, Phys. Rev. Lett. 70, 2845 (1993). 

[2] M. Gell-Mann, P. Ramond and R. Slansky, in Supergravity, edited by P. van Nieuwenhuizen 
and D. Freedman (North- Holland, Amsterdam, 1979); 

T. Yanagida, in Proceedings of the Workshop on Unified Theories and Baryon Number in 
the Universe, edited by O. Sawada and A. Sugamoto (KEK, 1979). 

[3] J. C. Pati and A. Salam, Phys. Rev. D 10, 235 (1974). 

[4] S. Dimopoulos, L. J. Hall and S. Raby, Phys. Rev. D 47, R3697 (1993); 
K. S. Babu and Q. Shafi, Phys. Lett. B 294, 235 (1992); 

K. S. Babu and Q. Shafi, Bartol Research Institute report BA-93-09 (1993), unpublished. 
[5] G. G. Branco and L. Lavoura, Phys. Lett. B 208, 123 (1988). 

[6] For the standard treatment of B^-Bd mixing, see for instance G. G. Branco and L. Lavoura, 
Phys. Rev. D 38, 2295 (1988). 

[7] J. Gasser and H. Leutwyler, Phys. Rep. 87, 77 (1982). 

[8] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); 

S. P. Mikheyev and A. Yu. Smirnov, Yad. Fiz. 42, 1441 (1985) [Sov. J. Nucl. Phys. 42, 
913 (1985)]. 

For a recent combined fit of the experimental data, see S. A. Bludman, N. Hata, D. G. 
Kennedy and P. G. Langacker, Phys. Rev. D 47, 2220 (1993). 

[9] N. Ushida et al, Phys. Rev. Lett. 57, 2897 (1986). 

[10] GHORUS GoUaboration, N. Armenise et al, GERN report GERN-SPSLG/90-42 (1990), 
unpublished; 



6 



NOMAD Collaboration, P. Astier et ai, CERN report CERN-SPSLC/91-21 (1991), un- 
published; 

Fermilab P803 and P860 proposals, unpublished. 

[11] V. Barger, R. J. N. Phillips and K. Whisnant, Phys. Rev. Lett. 65, 3084 (1990). 

[12] See B. Blumenfeld et ai, talk given at the Neutrino '90 Conference, CERN, 10-15 June 
1990, unpublished; and references therein. 

[13] T. K. Kuo and J. Pantaleone, Phys. Rev. Lett. 57, 1805 (1986), and Phys. Rev. D 35, 
3432 (1987). 

[14] The suggestion of this kind of double MSW effects was first made by L. Wolfenstein, Phys. 
Rev. D 45, R4365 (1992). 



7 



